Where is the charm quark energy loss at RHIC? 
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Heavy quark energy loss in a hot QCD plasma is computed taking into account the 
competing effects due to suppression of zeroth order gluon radiation below the plasma 
Q"^ I frequency and the enhancement of gluon radiation due to first order medium induced 

' Bremsstrahlung. The results suggest a surprising degree of cancellation between the two 

. medium effects for charm quarks and provides a possible explanation for the null effect 

' observed by PHENIX in the prompt electron spectrum in Au + Au at 130 AGeV. 
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Introduction 

^ ! Recent observations [I]-|H1 of large suppression moderate p_L ~ 5 GeV hadrons produced in 
pH \ Au + Au at the Relativistic Heavy Ion Collider (RHIC) have been interpreted as evidence for 
jet quenching of light quark and gluon jets [Oj-IIEl- Jet quenching was predicted [d HHj to 
occur due to medium induced radiative energy loss of high energy partons propagating through 
ultra-dense QCD matter. The quenching pattern therefore provides a novel tomographic tool 
that can be used to map out the evolution of the plasma density. Of course, even without final 
state interactions, gluon radiation associated with hard QCD processes softens considerably 
the lowest order jet spectra. This is taken into account through the (DGLAP) evolution of 
the hadronic fragmentation functions. Medium induced radiation is the extra gluon radiation 
that arises from higher twist final state interactions and depends on the optical thickness or 
opacity L/\ oi the medium. 

A similar quenching pattern was predicted fl9\ l20| l2T] to occur for heavy quark (c or 
b) jets. However, in [22j it was pointed out that the heavy quark mass leads to a "dead 
cone" effect for 9 < M/ E that reduces the induced radiative energy loss of heavy quarks as 
compared to light partons. Numerical estimates indicated that the quenching of charm quarks 
may be approximately about a half that of light quarks. Experimentally, PHENIX data [23] on 
"prompt" single electron production in Au + Au collisions at ^/s = 130 AGeV have provided a 
first test of heavy quark energy loss. Remarkably, no indication for a QCD medium effect was 
found within the admittedly large experimental errors (see also ref. In future runs, much 

higher statistics and wider pj- range will become accessible. 

In this letter we investigate whether the apparent null effect observed for heavy quark energy 
loss via single electrons could be due to other medium effects that compete with the induced 
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radiation. Our title is motivated by a similar puzzling null effect seen previously with light 
quarks at much lower SPS energies 

We study here the non-abelian analog of the Ter-Mikayelian [22j-j2H] effect. In j22j the 
suppression of radiation below a plasma frequency cutoff was estimated to be only ~ 10% 
effect on the induced energy loss. However, the Ter-Mikayelian effect on the zeroth order in 
opacity (L/X)^ radiation has not yet been considered quantitatively. The first estimates of the 
influence of a plasma frequency cutoff in QCD plasmas were reported in ref. [21] using a constant 
plasmon mass ljq [Sn]-|S2] ■ The k dependence of the gluon self energies and the magnitude of 
longitudinal radiation were not investigated. In this letter, we extend those results by taking 
both longitudinal as well as transverse modes consistently into account via the frequency and 
wavenumber dependent hard thermal (1-loop HTL) self energy Il'^'^{uj, k) P^-|37j. 

The dielectric properties of an isotropic plasma lead to a transverse gluon self energy 
IlTiyJ^k) with nT(ti;p;(0), 0) = oj'^iiQ) ~ /^^/3, where /x ~ gT is the Debye screening mass 
of a plasma at temperature T in the deconfined phase. In addition, long wavelength collective 
longitudinal gluon modes arise with ni(a;p;(0), 0) = ^^^^(O). This dynamical gluon mass sup- 
presses the radiation of soft u < ojpiik) gluons and shields the coUinear /c^ — > singularities 
that arise for massless quarks. 

The second part of our study is to generalize the GLV opacity expansion method [TT] 
to compute the first order induced energy loss including the kinematic effects due to both 
heavy quarks and massive transverse gluons. These mass effects change the formation times 
of the gluons and hence the destructive (LPM) interference pattern that reduces the radiative 
energy loss in comparison to the incoherent Bethe-Heitler limit. For ultra-relativistic light 
partons, the formation times are generally long compared to nuclear dimensions. The energy 
loss AE"*"*^ oc L? for very high energy light partons depends quadratically on the thickness of 
a static QCD medium j21 ^1 ^] , while for heavy quarks we show below that it is closer to a 
linear oc dependence. 

In this letter, we answer the question in the title by showing that there is a surprising 
degree of cancellation between the reduction of the zeroth order energy loss due to the QCD 
Ter-Mikayelian effect and the induced radiative GLV energy loss for heavy quarks. The technical 
details of the derivations and further discussion will be presented elsewhere |3H1 EH] • Our main 
result is that while the dead cone effect |22] is not sufficient by itself to explain the null effect 
by PHENIX [21], the approximate cancellation of medium effects for heavy quarks could be. 

Plasmon Effect on 0^^ order Energy Loss 

Our first goal was to test the accuracy of the commonly used simplifying assumption of neglect- 
ing the k dependence of effective plasmon mass as in jHl EHl IHI and also neglecting longitudinal 
modes. As shown below and in more detail in [SHj, this assumption was found to be accurate 
to ~ 10% as long as the asymptotic plasmon mass [SEj, i^oo = /^/ \/2 rather than the long wave- 
length ujpiik = 0) = ^2/3 uJoo is employed. The accuracy of the approximation also improves 
dramatically as the mass of the heavy quark increases. 

We computed the plasmon effect on the (zeroth order in opacity) vacuum energy energy 
loss, AE^J^^ = Jdk{dlT/dk + dlL/dk) using the optical theorem in the soft gluon limit ignoring 
the spin of the heavy quark. The zeroth order transverse and longitudinal radiated energy 
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losses per wave number for an initial quark jet with four momentum P = (-Ejp) are found to 
be given by 
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where = 2Pk + k"^ + is the invariant mass of the intermediate state, k"^ = u'^ — k^ and 
Pk = Euj — |p||k| cos 6' with uo = uxik) for the case of transverse radiation and uj = uji{k) 
for longitudinal radiation. The well known dispersion relations, urik) and ujiik) are obtained 
from u'^eriuj, /c) — k^ = and eii^j, k) = respectively. Analytic expressions for er and in 
the HTL or equivalently, semi-classical transport linear response approximations can be found 
in refs. [5^ IHHIITF] . The angular integration can be performed analytically if does not run. 
We perform the integration with running coupling at the scale as in conventional DGLAP 
evolution. 
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FIG 1. The reduction of the zeroth order (vacuum) energy loss for charm quark due 
to the QCD Ter-Mikayelian plasmon effect is shown as a function of the charm quark 
energy. The upper curve shows the vacuum energy loss if gluons are treated as massless 
and transversely polarized. The lower solid curve shows medium modified (but zeroth 
order in opacity) transverse fractional energy loss. The dashed curve shows the negligible 
additional effect of longitudinal plasmons. 



The numerical results for a medium characterized by a Debye screening scale /i = 0.5 GeV 
in Fig.l show that the longitudinal plasmon contribution to the energy loss is indeed negligible 
for the range of energies considered. However, the transverse plasmon mass effect reduces the 
zeroth order energy loss by, ~ 30% , relative to the vacuum case. The plasmon mass cutoff 
would effectively enhance the yield of high transverse momentum charm quarks were it not for 
the extra medium induced radiation discussed in the next section. 
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Medium induced radiative energy loss 



The second effect of the medium on the the propagation of heavy quarks is induced radiation 
caused by the multiple interactions with partons in the medium. In order to compute medium 
induced energy loss for heavy quarks we extended the GLV method [lU] to include the kinematic 
modification of the propagators. The new feature here is that unlike in previous applications, 
both the quark and and gluon plasmon mass (M and nig respectively) are taken into account. 

Since we found that the dielectric plasmon effect can be well approximated numerically 
neglecting the longitudinal modes and using the wavenumber independent asymptotic mass 
yu/v^ as the effective gluon mass, we will use this simplifying assumption in computing the 
induced radiative energy loss. 

The heavy quark generalization of the GLV energy loss to first order in the opacity L/X = 
J dxap is found to be [5^ : 



,(1) _ CpOis L , f fJ' dqj f dkj 9{xE - k±) 

2 L 2 , nj2^2 , _2w;^ 2 , n^2^2 , ^2\ 



2 ^^2^2 ^2 I (g^' - + + S) + + K (2) 

^ ' ^{{k^ - q^y + M2x2 + m2)((fcx + q±y + M^x^ + m2) ' 

where q± is the magnitude of the transverse momentum transfer between a target parton and a 
jet, k± is the magnitude of the transverse momentum of the radiated gluon, and A is the mean 
free path of the gluon. The simple analytic form of the destructive interference factor involving 
[AEx/LY arises for an assumed exponential distribution exp{—Az/L)/L between scattering 
centers. The expression in {■ ■ •} reduces in the M = rrig = limit to 2k'j_6{q± — kj_), and 
the GLV asymptotic result eq.(127) of [IHI can be recovered. We ignore the finite kinematic 
bounds on q± < ^/QET. The integral can then be performed analytically but is cumbersome. 
Qualitatively, the effect of increasing M is to reduce the relevance of the formation time factor 
(4£'x/L)2 in the denominator of the integrand. Formally, setting this factor to zero recovers 
the Bethe-Heitler result with AE^^^ oc OgEL/X modulo a logarithmic factor. 

The numerical results for the first order induced radiative energy loss are shown on Fig. 2 
for charm and bottom quarks. We fix the effective static plasma opacity to be L/A = 4 from 
the analysis of light quark quenching in Au + Au at 130 GeV from the results of Levai et al |41j . 
We assume here that a.^ = 0.3, /i = 0.5 GeV, and A = 1 fm for the plasma parameters. In 
the energy range, ~ 5 — 15 GeV, the Ter-Mikayelian effect reduces the induced energy loss 
in this extension of the GLV approach somewhat more than in the BDMS approximation [221 • 
However, on an absolute scale, this only corresponds to a change of S{AE^^y E) ~ 0.05, which 
is negligible. Note that with both dead cone and plasmon mass reduction, there remains a 
sizeable induced energy loss fraction AE^^^ / E ^ 0.2 for charm quarks while only about half 
that is predicted for bottom. Therefore, the dead cone effect by itself is not sufficient to explain 
the PHENIX heavy quark null effect. 
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FIG 2. The 1** order in opacity fractional energy loss for heavy quarks with (solid curves) 
and without (dashed curves) Ter-Mikayelian effect approximated by a constant rUg = fj,/V^ 
for transverse modes only. Upper curves correspond to charm, and lower to bottom quarks 
as a function of their energy in a plasma characterized by Og = 0.3, /x = 0.5 GeV, and 
L = 4A = 4 fm. 



The net heavy quark energy loss 

Fig. 3 shows the competition between the medium dependence of the induced energy loss and 
the zeroth order energy loss taking into account the Ter-Mikayelian effect. Even in the absence 
of a medium (L = 0), a charm quark with energy ~ 10 GeV suffers an average energy loss, 
AE^^^/ E ^ 1/3, due to the sudden change of the color current when it is formed in the vacuum. 
The dielectric plasmon effect reduces this to about AE^l^/E ^ 1/4. This contribution is 
independent of the thickness of the plasma as long as L is not too small. For very small 
L < l/rrig, the plasmon dispersion is smeared out due to the uncertainty principle, and AE^^l^ 
must approach AE^'^^ from below. 

The induced contribution, AE^^^, on the other hand, increases with L. Note that for charm 
quarks the thickness dependence is closer to the linear Bethe-Heitler like form, L^, than the 
asymptotic energy quadratic form^UlE]- More importantly, the two competing effects tend 
to cancel to a large extend for an effective thickness in the physically relevant range L. 
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FIG 3. The fractional energy loss for a 10 GeV charm quark is plotted versus the effective 
static thickness L of a plasma characterized by /x = 0.5 GeV and A = 1 fm. The dashed 
middle horizontal line corresponds to the energy loss in the vacuum taking into account 
the kinematic dead cone of radiation for heavy quarks j^. The lower horizontal solid 
line shows our estimate of the reduction of the zeroth order energy loss due to the QCD 
analog of the Ter-Mikayelian effect. The solid curve corresponds to the net energy loss, 

The net energy loss, AE^^l^ + AE^^^ is found to be remarkably close to the naive vacuum 
value AE^I^^ for the effective static opacity medium with L 3 fm and /i ~ 0.5 GeV, A ~ Ifm). 
These results suggest that in addition to the heavy quark dead cone effect, a high degree of 
cancellation between the competing medium effects may be an important factor that could 
explain the null effect found experimentally by PHENIXj2S] for Au + Au at 130 GeV. The 
different dependence of these effects on the plasma thickness L and on the transport properties, 
fi and A, should make it possible to test this explanation by varying the beam energy. A, and 
centrality at RHIC and LHC energies. 

In closing, we note a third mechanism that could also play a role in explaining the null 
effect reported by PHENIX. As shown in j32] , for quark momenta below ~ 3T detailed balance 
in a thermal medium gives rise to a positive feedback due to thermal gluon absorption that 
limits the energy loss. For a rapidly expanding hydrodynamic medium with transverse flow 
boost rapidity, rjT, even large energy loss cannot quench the spectrum below the blue shifted 
Pt ~ 3Te^'^ thermal distribution if the heavy quark comes into local equilibrium. As noted in 

, for charm quarks there may in fact be only a small room between the unquenched pQCD 
spectrum and the fully quenched hydrodynamical spectrum in the pt range accessible thusfar 
at RHIC. A decisive test of this explanation would be the observation of charm elliptic flow at 
RHIC. If, on the other hand, no charm elliptic flow is observed and the null effect persists to 
higher px, then the cancellation mechanism discussed above could solve the charm puzzle. 
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